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Here we write down and prove closed commutation rules for Fermi bilinears in a two- 
component Fermi system in terms of the relevant sea-bosons. We show how the commutation 
rules come out correctly within the RPA-approximation as do the zero temperature corre- 
lation functions. We write down hamiltonains of nonrelativistic interacting electron-hole 
systems and point out several attractive features such as the natural roles excitons play 
in this language. We then use this language to derive a set of equations analogous to the 
Semiconductor Bloch Equations in the presence of phonons and with external fields such 
as those present in pump-probe experiments. We solve these equations using parameters of 
interesting and topical materials such as GaN and compare with some recent experimental 
data. 



Here we write down the formal correspondence between Fermi bilinears and sea-bosons suitably generalised 
to two-component systems. The correspondence presented here is necessarily approximate but one that 
conforms to the spirit of the random-phase approximation, which in the two-component system corresponds 
to the exciton approximation. However, just as in the one-component case, new physics could be extracted 
by suitably generalising the notion of the random-phase approximation, we find that here too we may extract 
new physics by generalising the notion of the exciton approximation so that one goes beyond the dilute limit 
and includes experimentally relevent situations such as those in which a large number of real carriers are 
created by external fields. The discussion presented here follows closely the discussion in the response to 
the comment on our previous workEI. Let us write the Fermi bilinear sea-boson correspondence with spin. 



I. BILINEAR SEA-BOSON CORRESPONDENCE 
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Let us make the following identifications, 

CkT = Ck (3) 

Cki = dtk (4) 
Taking a cue from the one-component case let us now argue that( for both q = and q 7^ 0) 

c?-k-q/2Ck-q/2 ~ ^(1 - ne(k - q/2))(l - n,i(-k - q/2))aki(-q t) 

+ Vrie(k-q/2)n„(-k-q/2)aj^^ (q i) (5) 



Ai (k, q) = v/(l - rle(k - q/2))(l - n„(-k - q/2)) (6) 



A2 (k, q) = v/ne(k-q/2)n,,(-k-q/2) (7) 
The above form in Eq.( ^ automatically satisfies, 

[c'-k-q/2Ck-q/2,rf--k'-q72Ck'-q72] =0 (8) 

Next we would like the following to happen (for both q = and q 7^ 0), as it does in the Fermi language, 

[d-k-q/2Ck^q/2, = 
\,k'^q,q' - '^Lk'_q'/2^-k-q/2^k-q/2,k'-q72 " 4' -q72'^k-q/2^k+q/2,k' +q72 

- -^k.k'^q.q'll - (f^lk-q/2^-k-q/2) " (4-q/2Ck-q/2) ) (9) 

(1 - ne(k - q/2))(l - n,,(-k ~ q/2)) - ^^(k - q/2)n,,(-k - q/2) = f - (rf^_k-q/2'^-k-q/2) - (4-q/2Ck-q/2) 

(10) 

The expectation value is with respect to the full interacting ground state including (especially including) the 
external fields that allow for significant real populations to be generated. Let us now argue(again inspired 
by the one-component system). 



d^kC^-k = 4°^(k) - «k-q72^i(qi i)ak-q72^i(qi i) + X! "k+q72i (^H '^l )«k+qi /2i (QiCTi ) (11) 



qi(Tl qiCTi 

and, 



4ck = n'^eH^) + "k-qi/2ai(qi T)ak-q72^i(qi T) - XI «k+qi/2T (^1^1 )°^k+qi /2T (^1 ^1 ) (^2) 
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Here -^(k) and n), '{k) account for possible doping in the system. In other words, the presence of excess 
charge. These two put together obey the attractive identity (charge conservation), 

E - E ^-k'^-k = Q = E["e°Hk) - n^lk)] (13) 

k k k 

Now we move on to the off-diagonal(in the indices) parts(q ^ 0), 

^-k+q/2'^-k-q/2 = -A/»(k,q)aki(-q i) - A/«(k, -q)^^!^ J-) (1^) 

4+q/2Ck-q/2 = Ae(k, q)akT (-q T) + Ae(k, -q)4^(q t) (15) 

Assuming that the ground states (of the non-interacting system) are annhilated by the sea-bosons aki(q i) 
and akT(q T) 

aki(q i)\Free) = 0, akT(q ])\Free) = (16) 

This means, 

(4+q/2Ck-q/24_q/2Ck+q/2> = + q/2)(l " ^e(k " q/2)) = Al{k, q) (17) 
and similarly for the holes. 

Ae (k, q) = V«e(k + q/2)(l-ne(k-q/2)) (18) 

Ah{K q) = Vnfc(-k + q/2)(l - nhi-k - q/2)) (19) 

Then wc make a leap of faith and suggest that the same should hold even when there are inteactions present 
and even when external fields are present. That is, the f7,;i(k) — ((i^jj.(i_k) now represents the expectation 
value with respect to the full interacting ground state. Lastly we would like to point out the internal 
self-consistency of this approach by computing the commutator bettween the diagonal and the off-diagonal 
bilinears. We find much to our relief that no matter what the choices for the cofficients Ai, A2, Ae and 



are, we recover the following exact identities. 

[c?_k_q/2Ck-q/2>cJ,'^p] = '^-k-q/2Ck-q/2<^p,k-q/2 (20) 

M-k-q/2Ck-q/2,C^lp<^-p] = <^-k-q/2Ck-q/2<5p,k-|-q/2 (21) 

[4-Kq/2Ck-q/2, cj^Cp] = 4+q/2Ck-q/2 ('5p,k-q/2 " <5p,k-Hq/2) (22) 

[c?Lk+q/2'^-k-q/2) C^-p<^-p] = '^Lk+q/2'^-k-q/2 (^p,k+q/2 " ^p,k-q/2) (23) 
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Thus we have written down a potentially useful set of identities. Let us write down the hamiltonian of free 
electrons and holes. In the Fermi language it is, 

Hfree = e'(k)4ck + ^ e^k)(il^d_k + ^ fiiO&^^q (24) 
k k q 

The kinetic energy of the LO-phonon modes is also included. In the sea-boson language it may be expressed 
as follows : 

Hfree = " ^ (e^(k - q/2) + e'^Ck + q/2))4^(q i)akT(q i) 

kq^^O 

- E(£ + ^9)«kT(0 i)«kT(0 i) + ^ (e^k - q/2) - e^k + q/2))4^(q i)aki(q i) 

k ^ kq^tO 

+ ^ (e'^(k - q/2) + e^(k + q/2))4^(q T)aki(q T) + + Eg)ai^{0 t)aki(0 T) 

+ ^(e^(k + q/2)-e^(k-q/2))4^(qT)akT(qT) (25) 

kq^^O 

Again it may be noted that objects such as akt(0 t) and akj,(0 J,) are omitted from the formalism. The fact 
that the Fermi bilinears all evolve properly with respect to this hamiltonian is apparent without the need 
to perform any calculations. This is a strong indication that we are on the right track. It may puzzle the 
reader that we have included an object such as ak,j,(0, t) in the above formula. This is due to the following 
reason. The commutation rule [d_kCk, cj^''^^,^/] does not come out right if we don't. Let us now write down 
some typical interaction terms. 

He-h = -Y. 4+q/2'^-k'-q/2'^-k'+q/2Ck-q/2 (26) 

q#o k,k' 
This may be recast in the sea-boson language as follows, 

H,_n = ^^I^[Ai(k/2+k'/2+q/2,k'-k)a+^^^^,^^^^^^^(k-k' T)+A2(k/2+k'/2+q/2, k'-k)ak/2+k'/2+q/2T(k'-k i)] 

■37^0 k,k' 

[Ai (k/2 + k' /2 - q/2, k' - k)a^/2+k' /2-q/2i (k - k' T) + (k/2 + k' /2 - q/2, k' - k)a;^^^^^, (k' - k i)] 

(27) 

Let us now try and write down the e-e/h-h repulsion terms(let us now focus on an undoped system), 



He-e = Y.'^P^W^-^) (28) 



q#0 



if/.-^ = E^'°''(q)'^'(-q) (29) 

q#0 



^f'' = E + ^-q)(/''(q) - p'i^)) (32) 



p%q) = J2 Ae(k, q)akT (-q T) + E ^^C*' -q)4T T) (30) 

k k 

/^''(q) = - 1^ A^(k,q)ak^(-q i) - ^ A^(k, -q)4^(q i) (31) 

k k 

The coupling to phonons may be written as follows, 

It may be seen that only in the presence of real charge distributions do the electron-electron/hole- hole 
repuslion and coupling to phonons contribute appreciably to the hamiltonian. This means that in the 
undoped case in the absence of external fields we expect only the excitonic contribution broadened perhaps 
only via coupling to photons (which is ignored here). External fields, especially pump fields above the band 
gap cause significant real populations of carriers and these in turn relax by emitting phonons and through 
Coulomb scattering. Thus the formalism we have written down is simple and ideal for the study of these 
systems. The coupling to external fields may be written as, 

Hextit) = (^)le.t(i).p;c^[Ai(k,0)aki(0 T) + A2(k, 0)aj,^(0 i)] 



+ (U)l*^^(f).p-^^[Ai(k, 0)4^(0 T) + A2(k,0)akT(0 i)] (33) 

Let us now write down the various equations of motion of this system. 

i^aki(0 T) = (e^(k) + e''(k))aki(0 T) - ^ ^^^au-Qi(0 T) 

+ ^^^17^ [(1 - Ai(k,0)Ai(k- Q,0))ak_Qi(0 T) - Ai(k, 0)A2(k - Q, 0)at _qt(0 i)] + (^)^L«(t).p.cAi(k, 0) 

Q#0 

(34) 
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iJakT(0 i) = -(e=(k) + e'^(k))akT(0 i) 
- Yl ^^^^2 (k, 0) [Ai (k + Q, 0)4+Q^ (0 T) + A2 (k + Q, 0)ak+QT (0 i)] 



+ ( j^)Xxt W.p.cA2(k, 0) (35) 

i|auT(q T) = ^«kT(q T) + ^Ae(k, -q)p(^>(-q) + ^^qAe(k, -q) (36) 

i|aui(q i) = -^«ki(q i) - ^A,(k, -q)p('')(-q) + ^X,A,(k, -q) (37) 

»|4T(-qT) = ^4T(-qT) - ^Ae(k,q)p(^)(-q) - ^X,Ae(k,q) (38) 

= + ^A,(k,q)pW(-q) - ^X,A,(k,q) (39) 

i^Xq = (2iOz,o)P-q (40) 

4^-. = - i^(p^^'(-q) - /"H-q)) (41) 



The last four equations of motion only affect the electron and hole populations but do not impact directly 
upon the polarization or induced currents. This means that electron-electron and hole-hole repulsion and 
electron-phonon interaction change the distributions of electrons and holes and the electron-hole attraction 
determines the absorption spectrum. Let us now write down the electron and hole populations, 

n^(k) = -(at^(0 i)akT(0 i)) + (4^(0 T)aki(0 T)) 

" H^4-q/2i(q i)ak-q/2i(q i)) + Xl^4+q/2i(^ lW+q/2l{(l i)) (42) 

^e(k) = {al^iO T)aki(0 T)) - (4^(0 l)ay,^{0 i)) 
+ Z!^4-q/2T(q T)ak-q/2T(q T)) - X](4-Fq/2T(q T)ak-Kq/2T(q T)) (43) 

q/0 q#0 

In the above sets of equations we have ignored the contribution from objects such as akT(q i) with q 7^ 0. 
The reason being that these contributions are difficult to deal with. The practical consequences of this 
assumption means that we have to restric our attention to large k. Namely that we must ensure that k 
in the above equation is much larger than any inverse length-scale in the problem. Thus we expect our 
theory to be poor for k small. This is in fact the case as we shall soon find out. The analysis including 
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this large q effect Let us now introduce several propagators. These are going to be useful in ascertaining 
the influence carrier-carrier repulsion and carrier-phonon interactions have on the populations of electrons 
and holes. We can see from Eq.( ^2|) and Eq.( the manner in which the populations of electrons and 
holes evolve. If we ignore the terms that involve the sea-bosons akCT(qcr) thern we see that electrons and 
holes have the same distribution determined solely by the external fields. That is, so long as relaxation 
processess are ignored we find simple and intuitively appealing formulas for the distributions. It is also 
worth pointing out that this situation is entirely analogous(perhaps even equivalent) to the Semiconductor 
Bloch Equations(SBE)i. There we find that if one ignores terms beyond the Hartree-Fock approximation, 
then the populations of electrons and holes are indentical even though the effective masses are different. 
However once relaxation processes begin, the distributions respond appropriately and the system must be 
solved self-consistently. A quantity such as (aj^^(q(T)akCT(qcr)) is zero when the electrons and holes have 
"ideal" momentum distributions (that is, identically zero for undoped systems). However, when they start 
to acquire non-zero values due to external fields the above quantity also begins to acquire a non-zero value 
and the whole system proceeds to evolve accordingly. Let us now introduce several Green functions. 

Gn(k,k';q) = -i(TakT(qT,i)<.^(qT,0)) (44) 

Gi2(k,k';q) = -z(TaJ,^(-qt,i)4'^(qT,0)) (45) 

G2i(k,k';q) = ^t{T al^i-qi,t)al,^{q i,0)) (46) 

G22(k, k' ; q) = -z{T aki(q i, i)<.^(q i, 0)) (47) 

Gx.(k; q) = -i{X^it)al^iqa, 0)) (48) 

Gp,(k; q) = -i{P^^{t)al^{qa, 0)) (49) 
z^Gn (k, k' ; q, i) = S{t)5^ + ^Gn (k, k' ; q, i) - i^A,{k, -q) (T p(^) (-q, i)at (q t, 0)) 

+ ^Ae(k,-q)GxT(k';q,<) (50) 
V V 



i^G,2 (k, k' ; q, t) = ^Gi2(k, k' ; q, t) + z^t^A, (k, q)(r p(^) (-q, i)aj^, (q T, 0)) 
at nip V ^ I 



^Ae(k,q)GxT(k';q,i) (51) 
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z^G22(k,k';q,f) = <5(0(5k,k' - ^G22(k,k';q,t) +z^A;,(k,-q)(T p^^\-ci,t)al ^{q 1,0)) 

+ ^A^k,-q)Gxi(k';q,t) (52) 
i^G2i (k, k' ;q,t)= -^G2i (k, k' ; q, i) - i^T^ A^(k, q)(r pC^) (-q, i)a[, (q i, 0)) 

-^A^(k,q)Gxi(k';q,t) (53) 

i-Gxa(k;q,t) = {2inLo)Gp„{k;ci,t) (54) 

z|GpT(k;q,t) = (^)GxT(k;q,t) - ^(T p^-^H-q, i)a|,^(q T,0)) (55) 

i^Gp^{k;ci,t) = (^)Gxi(k;q,i) + ^(T pW(-q, i)a],^(q i,0)) (56) 

In order to solve this system we have to expand the Green funetions in terms of Matsiibara frequencies. We 
introduce a temperature just for ease of doing calculations. In the end we shall go to the zero temperature 

limit as this is the regime when the interpretations are the cleanest (here z„ = 27rn//3 and /3 = l/ksT ). 

]iS)Gn(k,k';q,^„) = -L,?^ - i!^i^Ae(k, -q)(T p(^)(-q, ^„)a;^, (q T, 0)) (57) 

{iZn - ^)Gi2(k, k' ; q, zn) = i^^^^A,{k, q)(r p(^)(-q, z.)aj^,^(q T, 0)) (58) 

{izr, + ^)G22(k, k'; q, z^) = ^-^k.k' + i^^h{K -<l){T P^'^H-Q. ^n)^, , (q i, 0)) (59) 

{zzn + ^^)G2i(k, k'; q, zn) = -i^%^A;,(k, q) (T p(")(-q, z„)<' , (Q i, 0)) (60) 
nih V ^ 

2 M^^Lo 

w(q, Zn) = v{q) - o , o2 (61) 

- i{T p(^) (-q, 0„)4,^ (q T, 0)) = ^ Ae(p, -q)Gii (p, k' ; q, Zn) + ^ Ae(p, q)Gi2 (p, k' ; q, Zn) (62) 

p p 

- i{T (-q, ^„)aj^, (q i , 0)) = - ^ A;, (p, -q)G22 (p, k' ; q, ^„) - ^ A^ (p, q)G2i (p, k' ; q, z^) (63) 
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(85) 



where, 



where, 



Therefore, 



p me 



.Wrn ,V ^ 1 + ^(q^^-) \- nW(p + q/2)-nW(p-q/2) 

^(q,z^„) = 1 + — 7^2^ T- pTq (68) 

^ p mft 



/ t / *x / / 1 ^'(q,2n) «A2(k, -q) 

(4,(q T)a.T(q T)) = (-^)(^) E e(e)(q,..„)(,,„_^). 



/ t ^ n / n\ ^ 1 v{ci,Zn) iA|(k,-q) 

(4,(q i)a.,(q i)) = ( — )(-) ^ (q, + 



(70) 



rrih ' 

Then when we go to the zero-temperature limit we have to integrate over all n. In Fig.l we see the pole 
structure of the above equations. Call i Zn = i z. Then, 

<<T(q T)a.T(q T)) = I dz ^(e)(q^,,) (71) 

In Fig.l wc sec the pole structure of the above contour integral. Let us assume that k.q > then the pole 
z = —i k.q/jTie is in the lower half-plane. We have to close the contour in such a way that this pole is 
excluded from consideration, since if it were included wc would have a formula for (akt(q t)ik|(q T)) rather 
than (aj(.|(q t)'ikT(q ?))• Therefore we have to close the contuor in the upper half-plane(and C = CI). 
Now, if we coimt the number of poles in the integrand we find that first of all. the zeros of the dielcc;tric 
function that lie on the positive imaginary axis of the z-plane(how many zeros are there, is an important 
question which we shall address subsequently) contribute. Then it seems at first sight that even the poles 
of ?}(q. z) contribute. The poles of this function lie at ±i Q.lo- However, upon closer examination we find 
that this is not the case. The poles of ^(q, z) are also poles of e^'^^(q, iz) and the two cancel. Thus the only 
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poles that contribute are the zeros of the dielectric function that lie on the positive imaginary axis. How 
many such zeros are there ? If one counts only the collective modes then one arrives at the conclusion that 
there are only two, one corresponding to the plasmon(modified by phonons) and the other corresponding 
to phonons (modified by Coulomb interactions). There is another mode that is equally important, indeed it 
would be a serious mistake to ignore this contribution, namely the particle-hole mode. We have encountered 
this problem before!. In our earlier article we presented an argumentS that shows how one may incorporate 
the particle-hole mode. In retrospect it seems that the approach presented there is not a good one, although 
it serves well to illustrate the importance of the particle-hole mode. Here we shall take the point of view 
that all energies are allowed as zeros of the dielectric function(for each q) but each comes with a weight 
corresponding to the strength of the dynamical structure factor at that energy. Thus for small q we recover 
naturally the collective modes but for larger q we start summing the particle-hole modes as well. There is 
really is no rigorous justification for this point of view except that it is physically well-motivated. 

/ t ^ t\ ( f\\ ('^\(^\ f i)(q,z) A^(k,-q) 

KT(q TKT(q T)) = ^ .(e)(q,..) (72) 



The poles are z = i uj'f\q)^ uj''j'\q) > satisfies e^'^^ci^ uj'f^ ) — 0. 

Let us first evaluate this. Similarly, one may write for holes, 

/ t ^ i\ / iN\ /1n/1\V^ v{ci,i ujf''') A?(k,-q) , ^, 

""'"^'"''""'^'■"^'^ AL. (.f^.)-^)^ (74) 

Let us now evaluate these quantities more explicitly. 

d d d 

-Q^\z=iuji^{'\,iz) = -P{^.,^z)\z=iuJI-Q^\z=iu,,v{(\, z) - v{ci, z)\z^i^j — \z=iujjP{ci.,iz) 

Id d 

\z=iuji^\z=iojiV{ci, z) - v{q^,z)\:,^i^, — \^^i^jP{ci,iz) 



v{(\, z) az oz 



Since, 



'■LO 



{i(q, z) ^ v{ci) - ,'7"„2'^ (75) 

Z ~T~ \ L T 



-U_,.(q,z)=^^,_^,^^, (77) 
^1 ■ r(a:z]- i ^ ^^LO^i >(q, »^j) s2 y - nk+q/2 
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As we pointed our just a while ago, it is necessary that we interpret the sum over / in a special manner so 
that we are able to recover both the collective as well as the particle-hole modes. The way this is done is 
through the following identification, 

/•OO 

where the weight is the dynamical structure factor normalised to unity. 

WM= (80) 

Jo d(jj S{q,uj) 

The dynamical structure factor is defined to be the the dynamical density-density correlation function fourier- 
transformed divided by the total number of particles. Let us first write down, 



<r T,«)> H.^, E ^-^^Snk) 

Since Im{t) e [0, -/?], if Im{t) < then, 

„zt «Ae(k',-q) 1 



(81) 



(p(«)(-q,^)at,^(qT,0)) = (-L)(A)y^^ 



iz 



iz) 



me 

If Im{t) > then, 

«.»,«, T.0).<'>(-<„<)> = 't*!^^',.,;,,) 

me 

me 

Here C+(C_) is the semi-circle in the upper(lower) half plane. Let us now take the complex conjugate of 
the above equation. Im(t) > implies Im{t*) < 0. Therefore if Im{t*) < 0, 

(.<->(,.*-)v,«.r,o))=-(i,/;; 4.^- 

me 

Define, 



(84) 



p(^''^)(q,0) = ^ Ae(k',-q)at,^(q r,0) (85) 

k' 

p(«''')(q,0) = ^Ae(k',q)ak'T(-qT,0) (86) 
k' 

therefore. 



11 



p(^)(q,0)=p('=''^)(q,0)+p('=^'')(q,0) (87) 

Then we have Im{t) < 0, 

(,(^)(-q,.),(^)(q,0)) = (f ) Prfxe- ^ ^(^e(k'-q/2) . (k'+q/2)) 1 

27r 7_oo ix~^ e(«)(q,ia;) 

k me 

For Im{t*) > 0, 

(p(^)(q,0)p(^)(-q,n) = (f ) Pc^— \ (gg) 



Combining these two, 



(rp(^)(-q,t)p(^)(q,0)) = (f ) P..e- E ^^'^^'^''^^'^".'f' ^^^'^W 



me 



)(q,ia;) 



Define the Green functior 



then, 



2^ 7-00 e(^)(q,za 



l?(q,i) = -(rp('=)(-q,t)p(^)(q,0)) (91) 



2?,et(q,c^)^^ (%;^' ; (92) 



The corresponding spectral function is the dynamical structure factor 



iVe 5(q, UJ) = -2Im{VreM, ^)) (93) 
e(,^)(q,c.) = l-i;'-(q,c^)P;(q,u;)+i;^(q,c.)PJ(q,c.) (94) 



ef)(q,o.) = -«'-(q,c.)P:(q,..) -t;'(q,c.)P:(q,c.) (95) 

This procedure ensures that we correctly incorporate both the collective (for small q) and the particle-hole 
modes. There is an alternative approach that comes to mind. That is the method of exact diagonalisation. 
Consider the hamiltonian, 

^' ^Yl — 4T(q t)akT(q T) - X! — i)aki(q i) ^Lob\b^ 

k,q k,q q 



u(q) 



V [Ae(k, q)auT (-q T) + Ae(k, -q)aj^^ (q t)] [Ae (k , -q)ak' t (q T) + Ae(k' , q)^, , (-q T)] 



Qt^O k,k' 
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E ^ l][A„(k,q)aki(-q i) + A^(k, -q)4^(q i)] [A;,(k' , -q)ak'^(q i) + A/.(k',q)aj^,^(-q i)] 



q^^o k,k' 
M, 



-^(&q + ?'lq)[Ae(k, q)akT(-q T) + Ae(k, -q)4^(q T) + A^k,q)aki(-q i) + A^k, -q)4i(q i)] 

(96) 

In order to diagonalise this we proceed as follows. Let us postulate the existence of dressesd sea-bosons 
d/o-(q) such that, 

h' =Y, ^iM)d\M)diM) (97) 

/,q,(T 

Now for some notation, a = e^h (correspondingly a =T,i, furthermore, S{a) is such that 5(1) = +1 and 

s{l) = -i). 

ak^iqa) = ^[ak^(qcr), 4^(q)]d7^(q) - ^[ai,^iq_a), di^{-q)]d\^i-q) (98) 
I I 

\ = J2K4A<l)]diM) - E[6q,rf7.(-q)]4(-q) (99) 

I,a 7,(7 

The inverse relation is, 

k k 

+ [di,{q), - [di,{q), b.^jbl^ (100) 



Therefore, 

uJia{q)di„{q) =^S{a) — [di„{q),al^{q_a)]ak„{cicr) - S{a) ^ — ^/^(q), ak<T(-qo-)]a],^(-qa-) 



+ ^ E ^^'^ C*' ^^^'^ ('l) ' Oka (-qa)] + A<, (k, -q) [di^ (q) , al^ (qa)]] [A,, (k' , -q)ak' ^ (qa)+ A,, (k' , q)aj^, ^ (-qa)] 
k,k' 

+ ^([d/^(q),6_q] + [d/<,(q),6jj])^[A^(k,-q)aka(qCT) + A^(k,q)a|,^(-q(T)] 
+ ^ (6q + 6Lq) ^(A, (k, q) [di, (q) , Ok. (-qa)] + A, (k, -q) [d/, (q) , (qa)] (101) 
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M<1) - S{a)^)[djM),al^{(ia)] = 5(a)^^^^A,(k, -q)p(-) (q, /) (102) 

(a;/.(q) - 5(a)^^)[d/.(q), aka(-qa)] = -5(a)^%^A,(k, q)p(-) (q, 7) (103) 

p(«)(q, i) = Y. ^^(k' q)M/(q)> «kT(-q T)] + ^-C*' -q)['^/(q)' 4T(q t)] (io4) 

k k 

p('*)(q,/) = -^A^(k,q)[dj(q),aki(-qi)]-^A^(k,-q)[dz(q),ai,^(qi)] (105) 



(a;,,(q) - Oio)[rf7.(q), b^,] = 5(<T)^p(-)(q, 7) (106) 

V 1/ 



{u;iM) + ^Lo)[diM),b-^]=-S{<y)^P^'\ci,I) (107) 
Prom this we have the following fact that are zeros of the dielectric function, 

e(^)(q,4"^) =0 (108) 

e('')(q,a;f^) =0 (109) 

e(^'''Hq,w) = 1 - v{ci,uj)P^^'''\ci,uj) (110) 

and, 

^;(q,a;)=^;(q)+ (111) 
and P(^'^)(q, w) is the usual RPA-polarization bubble. If we now make use of the fact that 

4( 



[d/cr(q), rf|^(q)] = 1 then we have. 



El['^^-(q)'«L(q^)]|' - El[^^-(q)'«k<.(-q^)]l' 

k k 

+ |[d,<,(q),6t]|2 _ |[d,^(q),6_q]|2 = i (112) 
^ K (k, -q) - A^ (k, q) ^ ^(^;| ^,12 ^ ^(q. -^t7) ^2 



+ f [p^-Hq,/)] \jt^^"f,. =l (113) 



Therefore, 
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^ ^^'^ ^ V K.(q) + ^)^ ^ V {u:],s^)-niy ^^^^^ 



Therefore, 



/ y^I ,h vqj j 



After some algebra it is clear that Eqs.( [73|) and ( |7J) are identical to Eqs.( 117) and ( 118) respectively. 
Let us now solve the fundamental equations namely Eq. ( ^) and Eq. ( ^ ) . For this we first would like to 
decompose the various fields in the exciton basis. 

aki(OT) = E^^We— ^ (119) 

akT(Oi)=e*('^('^)+^'(''»*akT(Oi) (120) 
This means we may rewrite these equations as follows : 

+ E ^^^lE - AiC^' 0)^1 - Q' - Q)e'(^^-^'^) 'Dj{t) 

Q#0 k,J 

- E ^Kk)Ai(k, 0)A2(k - Q, 0)e"^ *e-^('^(''-^)+^'(''-'^»*4_Q^(0 i)] (121) 

k 

45kT(0 i) = - E ^^A2(k,0)[Ai(k + Q,0)(^}(k + Q)e^^.^ *,-.(.= (k)+.'-(k)) t^t 
Qt^o 

+ A2(k + Q,0)e-'(^^('^)+^"('^)) tg,(.= (k+Q)+.'-(k+Q)) (0 j)] + (\A)A,,^^).p,, A2(k,0)e-'(^^('')+^"('')) * 

/iC 

(122) 
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dt lie ^ 



E ""'"^^ - Ai(k,0)Ai(k',0))^o(k')Do(i) 



kVk 



(^5(k)Ai(k,0)A2(k',0)e^(^°-v) *a|^,^(0 i)] 



^^^k(i) = (^)e^^'' *^L,(i).p-;cAi(k,0) 



^ ^'^^('^^ k ) _ ^^^^^ o)A,(k', 0))yo(k>'(^--^°) *Do(i) 

kVk 



kVk 



-Ai(k,0)A2(k',0)e'(^--v)*a;^,^(0 i)] 
^^4t(0 i)= E '''''^V ^ ^ A2(k,0)Ai(k',0)yo(k>-^(-°---) *^oW 

kVk 



+ E ^^2^^^^^7^A2(k,0)[Ai(k',0)e-^(V--) 

kVk 



+ A2(k',0)e^(^''-V)*a|^,^(0 i)] - {}3.^)Al,,{t)-Pvc A2(k,0)e">' 



<^o(k) = ^(47r) 



1 2 /qx 

3 (^ In"^ _L J.2\2 



(k) = no(k) + y r'^\;'^' yne{\^ - q)(l - ne(k)) 



2 



E(^^)V(k)(l-fle(k + q))- (Z'^^^^-^'^))^ 
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no(k) = |¥'o(k)P(^iW^oW) + (^JcW^k(i)) 
+ M^)e'^'''-'°^'{Di{t)Do{t)) + M^)e-'^'''-'°^' {Dimwit)) - (a|,^(0 l)a^^{0 i)) 

no(k) = |<po(k)|2((^5(i)^5(i)) + {DmOhm + (^kW^kW) + (^UO^LW) 
+2^o(k)cos((ek - eo)i)[(^i;(i)^5W) + {DUtMm 
2<^o(k)sm((ek - eo)i)[(^i;(i)^SW) - (^UO^SW)] - («kT(0 i)«kT(0 i)) " (4t(0 i)4T(0 i)) 

ne(k) = no(k)Ae(k) + (1 - no(k))Se(k) 
n/,(k) = no(k)A(k) + (1 - no(k))Bft(k) 

^ i+Tf (k) 
Be(k) = ^ 



T,«(k) 



rf (k) = Y^CAr^fM^ - q) 

Ke(q) + le-i)' 

Ke(q) + ^ + 2^;)^ 



Afc(k) 



s^k) 



r,'>(k) 

^ l+T^(k) 



1 , ^^'^'.''(k) 
T^(k) 
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Th^) = EC-^^^ni - n,0^ - q)) (139) 



T.^(k)^E(^)V(k..)^p,^^ (140) 
+ E ^ ^ Mm^ - Ai(k,0)Ai(k ,0))^o(k )J^(t) 

kVk 

+ E "'^'^^^i ^^ smiieo - 61,0 t)[(^o(k)(l - Ai(k,0)Ai(k',0))J-,(t) 

kVk 

-^o(k)Ai(k,0)A2(k',0)a;;,^(Oi)] 
+ E ^-^(^^~^) coa((eo-ei,0 i)[^o(k)(l-Ai(k,0)Ai(k',0))D^,(i) 

+ ^o(k)Ai(k, 0)A2(k', 0)a^,^(0 i)] (141) 



l^^(t) = -(M)Ai(eo,t)p.cE^i(l^'0)^o(k) 



^ !^£^i(k_JO^^(k)(l-Ai(k,0)Ai(k',0))^o(k')^S(^) 



^ !Wk_Jilco5((,^_,^,) f)[^o(k)(l-Ai(k,0)Ai(k',0))^;;,(i) 

kVk 

-<^o(k)Ai(k,0)A2(k',0)a;^,^(Oi)] 



^ !Wk_kl^.^(^^_^^,) t)[<^o(k)(l- Ai(k,0)Ai(k',0))D^,(i) 

kVk 



+ ^o(k)Ai(k, 0)A2(k , 0)dl,^iO I)] (142) 



^^Dlit) = -(M)A^(ek,0P.cAi(k,O) 
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+ Yl ^^2!L(]i_Jil(i _ Ai(k, 0)Ai(k', 0))<^o(k')Nn((ek - eo) t)^5(i) + cos((ek - eo) tM{t)] 

-Ai(k,0)A2(k',0){sm((ek-ek') i)ak'|(0 i) - cos((ek - e^') i^'T^^ (1^^) 
-§^Di{t) = (M)^5,(ek,i)p.cAi(k,0) 

+ E - Ai(k,0)Ai(k',0))^o(k'){cos((ek - eo) t)^S(t) - sm((ek - eo) i)^^(i)} 

kVk 

+ E - Ai(k,0)Ai(k',0)){cos((ek - e^,) t)Dl,{t) - szn{{e^ - e^,) t)Dl,{t)} 

-Ai(k,0)A2(k',0){cos((ek-ek') t)al^{0 I) + sin{{ek - e^,) t)a^/^(0 i)] (144) 
^"LtCO i) = E '''"^V ^ ^ A2(k,0)Ai(k ,0)^o(k'){co5((eo - ek) t)^S(t) + sm((eo - ek) t)^^(t)} 

kVk 

+ E '''"^V ^ V 2(k,0)[Ai(k',0){co3((ek' - ek) t)Dl^,{t) + sin{{e^^ - ek) i)^t'W} 

kVk 

+ A2(k',0){cos((ek - ek')i)«k'T(0 + - e^')t)dl,^{0 i)] - (|^)A5f (ek, i)p.c A2(k,0) (145) 

^«LT(Oi) = E "'"^V ^ ^2(k,0)Ai(k',0)^o(k ){-szn((eo- ek) t)^5(t) + cos((eo - ek) t)^SW} 

kVk 

+ E ^^^^^^^;F^A2(k,0)[Ai(k',0){-sin((ek' - ek) t)^ (t) + cos{{e^> - e^) t)Dl,{t)} 

kVk 

+A2(k',0){sm((ek - ek')i)«k'T(0 ^) ' - ek')i)5L'^(0 i)}] 

+ (M)^V(ek,i)p.c A2(k,0) (146) 
In order to simplify the calculations further, let us define, 
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2m' ^ 



m 2m ^ 



m 2m 



2 m 



m 



5(q,c^) 



PKq,w 2) 



P.(q,w) 



(1 - z;(q, c.)P,.(q, c.))^ + i;2(q, c.)(F,(q, c.))^ 



(147) 
(148) 
(149) 



p(q,t^) = {-( 



4a;a; 



(150) 



II. COMPARISON WITH SEMICONDUCTOR BLOCK EQUATIONS 



The equations presented in the previous section namely Eqs.( |121| ) and ( |122| ) are intended as alternatives 
to the usual Semiconductor Bloch equations used to study semiconductors. Let us now try to solve this 
system with a pulse field. That is, we apply an external field with central frequency lox and assume it lasts 
for a time tx and therefore we have a spread in frequency Tx — ^tx Jtx- We would like to see the evolution 
of the polarization and populations in this case. This exercise also enables us to compare our results with 
those of the SBEs and ascertain where the differences lie. To this end let us set, 

X.t(i) =P.cAx(0,t) (151) 

where, 

Ax{E,t) = Ao r dc. ^-2ih^e<--^) * (152) 

J -00 (w - wx)^ + ^x 

For comparison the SBE is reproduced below. 

ghhi^t) = mft(k) = z(dLkd-k), = i f{k)gee{kt) = me(k) = i(4ck) 

ghei^t) = i(d_k(Ock(i)) = i p(k) 

i^PhhCkt) = 2 Re{n(kt)gl^{kt)) + i?w.(ki) 
ot 

i^Qheikt) = -n{kt){i - 2ghhikt)) + (ehCk) + e,(k) - 2S(ki))5,,e(ki) + Rhe(kt) 



df 



W = |^A:,,(t).pvc -*^t'k-k'5'^e(k't) 

k' 
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Define, 



E{kt) = -iJ2vu-u'9hhO<:t) (153) 

k' 



S/(k) 



dt 

dt 

dpijk) 

dt 



2^1fl(k,i)i)/(k) - 2^1/(k,t)pH(k) (155) 
= f2fl(k,i)(l - 2/(k)) - 2E{k,t)pi(k,t) (156) 
= i7/(k, t)(l - 2/(k)) + 2S(k, t)pR{k, t) (157) 



fjjj(k,i) = (^)iext(i)p^c cos{{k'^/2fx + Eg-ujx)t) + ^ w^.k' N(k')cos((A:2 - k'^)t/2n) - pi{k)sin{{k'^ - k"^)t/2fj,)] 

(158) 

ni{k,t) = ( — )iext(i)P^c sin{{e/2fx + Eg-Lvx)t) + ^ t;k-k' [P/(k')cos((A;^ - ^'2)^/2^*) +p/{(k')sin((fc2 - k'^)t/2ij)] 

kVk 

(159) 

S(k,i)= ^«k_k,/(k') (160) 

kVk 



III. OPTICAL CONDUCTIVITY 

First define the total polarization, 

P{t) = ^(d_kCk) (161) 

k 

From tliis we may obtain tiie Fourier component, 



P{uj) = / dt P{t) e«^* (162) 
J —00 

Now since, 

Ar{t) = Ar e-' '"^ ''5{t - t) (163) 

E{t) = - ^.e--- r^'(t-r)^p'^ e-^-'E{u) (164) 

21 



E{oj) = {ioj)Ar e*('"-'"^'^ (165) 

Since, 

jM=p„eJ^H (166) 

We have 

a{u)=pyce LtniA^^o tt^t^ (167) 
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Here, 

H^) = Y.f '^^ ^(1 - ^e(k))(l - nft(k))[^o(k)^^)(i)e*("-^°)* +4-)(0e*("-^'')*] 
+ Y.f 'ii V«e(k)n;,(k)e*('^-^'')*4^(0, i, -) 
V(l - ^e(k))(l - n,.(k))[</Po(k)4+\i)e'^"-^°^* + D^^\t)e'^^-'-^'] 

K 

+ ^ / Vne(k)n^(k)e*('^-^'')*4^(0, i, +) (168) 
k •''^ 

i{Dl+\T) - D^o~\r)) = (^)e'(^''-'^-)M.p.eE<^o(k)Ai(k,0) (169) 

i{Dl+\T) - Di-\T)) = (M)e^(i+^''-"-)M,p,eAi(k,0) (170) 
i(akT(0 i, +) - akT(0 i, -)) = (M)e-^(i+^''-'^-)M,p„,A2(k, 0) (171) 
P^(^)=l J^dl" V(l-rie(k))(l-n,(k)) 

[(po(k)£'^"^'^(Ocos((w-eo)0-¥'o(k)^^"^'^(t)sm((w-eo)0+£'k"^'^(Ocos((a;-ek)0-^k~^'^(Os*" 
+ j -^2^^^^ J dt ^/nJk)KiJk)cos{{uJ - ek)t)a^^{0,l,-) 

+ y (27r)3* J '^^ -\/rie(k)n/i(k)sm((w - e\^)t)al,^{0, [, -) 



+ X Wr'^^ Jr '^^ ^(l-"e(k))(l-n,(k)) 



[(po(k)D^+)'^(t)cos((w-eo)0-Vo(k)^^+^'^(Osm((a;-eo)0+^i+^''*(Ocos((w-ek)0-^k^^'^(*)s'»"^ 

r-oo 4^^2 ^oo 



(2^) 



2 (-oo 

jdk J dt ^/neO^)nh{^)cos{{LJ - €k)t)a^-^{0, 1, +) 



+ J J2^dk J dt y^nJXjnJy}sin{{io-e^)t)dl^iO,l,+) (172) 
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+ / TTT-^T* V^^k)Mk)sm((w - ek)t)akT(0'i.-) 

-/ TTT-^dk I rft -v/^^I(k)^!^(k)cos((a; - ek)t)ak|(0,i,-) 
Jo (.^""J J-oo 

+ io (2^'^^ X V(l - ^e(k))(l - n,(k)) 
[<^o(k)^^+^'''(t)sm((a;-eo)0+V'o(k)^^+^''(i)cos((a;-eo)i)+£>^+^'''(t)sm((a;-ek)i)+^^^ 
+ ^^rffc y V^e(k)n,,(k)sm((a; - ey^)t)a^^{Q, i, +) 

"io (2^'^'' J V^^O^)Mk)cos((a;-ek)i)akT(0.i,+) (173) 

^(+^'«(r) - ^^^'''(t) = (M)A,p^^si„((eo - a;x)r) H ^dk ^o(k)Ai(k, 0) (174) 

^(+)'^(t) - D^-^'^t) = -{\^)AMOs{{eo - u;x)t) H ^dk (^o(k)Ai(k, 0) (175) 

4+)'«(r) - 4-)'^(t) = (M)A.p,eSm((ek - a;x)r)Ai(k, 0) (176) 

Di+^'\T) - dI-^'\t) = -(M)A,p,,cos((ek - a;;,)T)Ai(k, 0) (177) 

akT(0, i, +) - akT(0, i, -) = -(|^)A^p„,sm((ek - a;x)r)A2(k, 0) (178) 

akT(0, i, +) - akT(0, i, -) = -(|^)A^p,;cCos((ek - a;x)T)A2(k, 0) (179) 

P{uj)=Y^ r dt V(l - ne(k))(l - nh{mM^)Di~\t)e'^^-'°^' + ^-^Oe'^""^"^*] 
k 

+ Y.f rft V^e(k)nft(k)e'('^-^'')*4^(0, i, -) 

+E / V(i - ne(k))(i - nh{mM^W\ty'^''-'"^' + 4+)(0e*("-^'')*] 

k -^^ 

Vne(k)%(k)e*(-^'')*4^(0, i, +) (180) 

K 
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The Optical Conductivity 

cos{{lj — u]x)t)Pi{uj, At) — sin{{uj — lux)t)Pr{lj, Ar) — cos{{lo — LOx)r)Pi{uj, 0) + sin{{LU — lox)t)Pr{lu, 0) 



Re{o{Lo)) = je| LtA^-^o- 



LuAr 

(181) 



A quick check of dimensions. In units of 7i = c = 1, all masses are of dimension inverse length we take to be 
centeimeter. All times are in centimeters. Charge is dimensionless. Since ^-^AtPvc t is dimensionless, it follows that 
[Ar] = L~^. It is easy to check that [Re{a{u!))] = L"^ as it should be. 



IV. RESULTS AND DISCUSSION 

The equations written down in the previous sections have to be solved numerically. It is worthwhile to 
point out some pitfalls and problems. Let us first focus on the SBEs. As we pointed out in our earlier work0 
involving the SBE's the sums over k have to be carried out in a special manner so as to avoid potential 
divergences. 

' (fc' + k) 



(k' - k) 

{J^max H~ ^) 



(182) 



Further, it was suggested by Binder et.al.i that we should use a momentum cutoff kmax = ^'2/ dBohr, where 
otBohr is the exciton Bohr radius. The justification for this stems from the fact that beyond this cutoff the 
probability of the electron existing is negligible. This assertion is true only if the pump field frequency is 
below or equal the band-gap. When the frequency is well above the band gap the situation is less clear and 
care must be taken inorder not to lose features that may be present at high momenta. Unfortunately we 
have found that even when the pump field has a frequency equal to the band gap the prescription of Binder 
et.ai has some problems. In particular, we have found that if we try and sneak a peek at the form of the 
distribution for fc >> 12/ aBohr^ we find a periodic pattern suggesting therefore that electrons can exist at 
(arbitrarily) high momenta long after a pump field whose frequency is at the band gap is switched off (after 
a time long enough so that we may still meaningfully talk of a well-defined frequency) . This is a paradoxical 
and counterintutive result that has been gloosed over by the pioneersi. The arbitrary cutoff of Binder et.al. 
shold not be taken too seriously. In order to make more sense out of all this we have to claim that the SBE's 
produce the correct momentum distributions only for small enough fc, and we have to use some judgement 
as to where we should cutoff the distributions. 

The sea-boson analogs of the SBE's written down above have their own numerical problems. First is the 
fact that even in the two-component case the sea-boson technique works well only when q << k. Since we 
have chosen to study only q = 0, parts of the hamiltonian, it seems that we are in good shape. However we 
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find that even then there is a cutoff small k = kmin cutoff below which the momentum distributions become 
unphysical (larger than unity). This is true if we use the formula in Eq.( |130| ) n^Ck) ~ «/i(k) = no(k). 
Further we find that this identification is the analog of the SBE. It is comforting to know that the sea-boson 
technique is equivalent to the SBE in some limit. In the SBEs, the momentum distributions of the electron 
and holes are identical even if the effective masses are very different. This is due to the fact that the SBEs 
neglect the collision terms responsible for the asymmetry that we would otherwise expect. Similarly, the 



sea-boson equations at the level of Eq.( 13C) neglect the repulsion and phonon terms. However the SBEs 
do include repulsion at the Hartree-Fock level, therefore the analogy between the two is not exact. In Fig.l 
we see how far we may take this analogy between the SBEs and the sea-boson equations. The approach 
toward unphysical behaviour for the small k limit of the momentum distribution obtained using the sea- 
boson equations is also seen. When we include the effects of repulsion and phonons, the answers change 
quite dramatically. In fact they are so very different from the SBE results that we have decided not to 
publish them. It will take some more time before a thorough analysis is completed and all the ramifications 
are explored. For now we shall assume that the momentum distribution is that given by no(k) or that given 
by the SBEs. Let us first write down some formulas that relate the real part of the conductivity to the 
absorption coefficient. We may expect the two to have qualitatively similar features. However just to be 
sure and so that we don't make any mistakes having gotten this far, let us write down the formulas. They 
are a combination of the formulas from the text by Haug and Kochi and the one by ManahH. The transverse 
dielectric function may be decomposed as follows. 

e{uj) = ei{uj) + i e2{uj) (183) 

The abosrption coefficient, refractive index and the real part of the conductivity are given by(in units 

a{uj) = -5^e2(w) (184) 
n(U!) 

n{co) = {i(ei(c.) + ^el{cj) + elicj))}i (185) 
The real part of the conductivity is, 

Re{<j{ujj) = ^e2{Lu) (186) 

It is better not to use the Kramers-Kronig relations as our answer for the real part of the conductivity is 
undetermined upto a factor (actually it has no reason to, it just so happens that the magnitude does not 
agree with observations). We may write down a formula for the imaginary part of the conductivity just as 
we did the real part. 

-sin{(i^ ~ i^x)t)Pi{i^,At) ~ cos{{uj - ujx)t)Pr{i^,At) + sin({Lu - ljx)t)Pi{uj,0) + cos{(i^ - i^x)t)P 
lm(a(uj)) — \e\LtA^^o 



(187) 
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/m(aM) = - — (eiM-1) (188) 

Therefore we may deduce the optical dielectric function and hence the absorption coefRcient. All this would 
not be neceesary if the refractive index was close to unity and then 62 is negligible in comparison with ei for 
all frequencies. Then the real part of the conductivity would be proportional to the absorption coefRcient. 
Let us compare the experimental magnitude of the absorption coefficient and the energy lo. We find according 
to the experiments of Song's group!, |a| ~ lif'cm^^ whereas lo ~ 2.07r/(352 x lO^^cm) = 1.78 x lO^cm^^. 
we can see that these two quantitcs are comparable to each other suggesting thereby that the real part of 
the dielectric function is not close to unity. 
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